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Abstract xThe effects of the symmetry potential and the isospin dependent 
. in-medium nucleon-nucleon (NN) cross section on the number of proton(neutron) 

emissions N p (N n ) are studied respectively within an isospin-dependent quantum 
molecular dynamics (IQMD) model. The isospin dependent in-medium NN cross section 
is found to have a strong influence on N p (N n ) but N p (N n ) is not sensitive to the 
symmetry potential for the neutron-deficient colliding system at relatively high 
energies. We propose to make use of the N p (N probe to extract information on the 

"q ■ isospin dependent in-medium NN cross section. 
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1 Introduction 

With the rapid advance in radioactive beam physics a better understanding of the 
isospin degree of freedom in nuclear collision dynamics may provide us with useful hints 

on how to extract reliable information about the in-medium NN cross section and the 
symmetry potential [1-5]. To obtain this information several interesting isospin effects in 

heavy ion collisions have been explored both experimentally [6-16] and theoretically, 
e.g., [17-25] over the last few years. Bao- An Li et al. [19,25] have found that information on 
the symmetry potential can be extracted by studying the neutron to proton ratio of 
preeqilibrium emissions in the relative low energy heavy ion collisions (HlC).In 

contrast, in this paper our studies within IQMD found that the numbers of 
proton (neutron) emissions N p (N n ) depend sensitively on the isospin dependent 
in-medium NN cross section and weakly on the symmetry potential in the relatively high 
beam energy region for the neutron-deficient systems. In this case,N p (N ri ) is a probe for 
extracting information about the isospin dependence of in-medium NN cross section in 
the relative high beam energy region for the neutron-deficient systems. 
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2 Theoretical Model 



The dynamics of intermediate energy heavy ion collisions described by QMD con- 
tains two ingredients: density dependent mean field and in-medium N-N cross section. To 
describe isospin effects appropriately, QMD should be modified properly: the density de- 
pendent mean field should contain the correct isospin terms including symmetry potential 
and Coulomb potential,the in-medium N-N cross section should be different for neutron- 
neutron ( proton-proton ) and neutron-proton collisions in which Pauli blocking should 
be counted by distinguishing neutrons from protons. In addition, the initial condition of 
the ground state of two colliding nuclei should also contain isospin information. The main 
physics ingredients and their numerical realization in the IQMD model can be found in 
Refs. [3, 28, 30, 31, 34]. In the IQMD model, the density distributions of colliding nuclei were 
from the calculations of the Skyrme- Hatree-Fock with parameter set SKM* ' 32 1 and the 
initial code of IQMD was used to determine the ground state properties of the colliding 
nuclei, such as the binding energies and rms radii which are the same as the experimental 
data.In the presence calculations the interaction potential in the IQMD were determined 
as follows: 

U{p) = U Sky + V Coul + U sym + V Yuk + JJ MDI + U Pauli (1) 
jjSky - g d ens ity-dependent Skyrme potential, 

U Sky = a{—) + /?(— ) 7 (2) 
Po Po 

V c is Coulomb potential. U Yuk is the Yukawa potential , 

U*<* = t 3 exp(^M)/^^ (3) 
m m 



jjMDI j g m omentum dependent interaction, 



U MDI = U ln 2 Mrt-m) 2 + l]- (4) 

Po 



jjPauli ig the p auli p otentia i j 



U PauU = Vvi^-fexA-^-^ - { ^-pl}5 p%Pj (5) 



PoQo 1% 2pj 



o 



I 1 for neutron-neutron or proton-proton 
ViV] j o f or neutron-proton 

jjsym - g ^ e symmetry potential. In the present calculation, two different density dependent 
symmetry potentials I 2,3 ' 29 ] are used, i.e., Ul ym = cF\(u)5t z and U^™ 1 = cF2(u)[5t z + \5 2 }, 
where t z = 1 for neutron and t z = — 1 for proton, Fi(u) = u and ^(n) = u 2 , u = p/po- 
S is the relative neutron excess 5 = ^ n ~^ p = Pn ~ pp .Here c is the strength of symmetry 

potential, taking the value of 32 or MeV (the c=0.0 case is denoted by Uq 1 ™). p,p , 
p n and p p are total density and its normal value, neutron density and proton density, 
respectively. It is worth mentioning that the recent studies on collective flow in HIC at 
intermediate energies have indicated a reduction of in-medium N-N cross sections. An 
empirical expression of the density dependent in-medium N-N cross section! 26 ! is given by 

= {1 + a±)ogZ (6) 
Po 
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Where the parameter a « —0.2 has been found to reproduce the flow data. <J^^ is the 
experimental N-N cross section' 27 ! . The parameters of the interaction potentials are given 
in table 1. 
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Table l.The parameters of the interaction potentials 

The free neutron-proton cross section cr^ e is about a factor of 3 times larger than the 
free proton-proton or free neutron-neutron cross section below about 400 MeV (in the 
Lab). It is worth to mentioning that the relationship between the neutron-proton cross 
section and neutron-neutron (proton-proton) cross section depends also on the modifica- 
tion of the nuclear density distributions during reactions. We construct clusters by means 
of the isospin-dependent modified coalescence model ' 34 1 ,in which particles with relative 
momentum smaller than po= 300MeV/c and relative distance smaller than Rq= 3.5 fm 
are used. We make use of the restructured aggregation model' 35 ' to avoid the nonphysical 
clusters after constructing the clusters, until there are not any nonphysical clusters to be 
produced. 

3 Results and Discussions 



3.1 Checking the IQMD model 

In order to check the IQMD code with the above parameters, the multiplicity of the 
intermediate mass fragments N,i m f for the reactions 58 Fe + 58 Fe and 5S Ni + 58 Ni at 
the beam energy E = 75MeV/u has been calculated by using the IQMD code with the 
above parameters. The multiplicity of the intermediate mass fragments (IMFs) is defined 
as the number of fragments with charge numbers from 3 to 18. The calculated results 
are compared with the experimental data ^ on the same scale in Fig.l which gives 
the correlation between the mean value of the intermediate mass fragment multiplicity 
Ni m f and the charged particle multiplicity N c . The solid (open) circles represent the 
experimental data for the reaction 58 Ni + 58 Ni ( 58 Fe + 58 Fe) at E = 75MeV/u and the 
solid line (dot line) denotes the IQMD results for 58 iVi + 58 Ni ( 58 Fe + 58 Fe).It is clear 
that the present IQMD predictions are in satisfactory agreement with general features of 
the experimental data which means that IQMD with the above parameters is a reasonable 
transfer theoretical model for simulating the dynamical process in intermediate energy 
heavy ion collisions. 

3.2 A probe of isospin dependent in-medium NN cross section by nu- 
clceon emissions 

The isospin effects of the in-medium NN cross section on the physical quantities arise from 
the difference between isospin dependent in-medium NN cross section denoted by cr lso in 
which a np > a nn =a pp and isospin independent NN cross section denoted by a nmso in which 
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&np = &nn = &pp- Here a np , a nn and a pp are the neutron-proton ,neutron- neutron and 
proton-proton cross sections respectively. 

Here N p (N n ) includes all of protons (neutrons) emitted during the nuclear reaction. To 
identify free nucleons,a phase-space coalescence method has been used at 200fm/c (when 
N P (N n ) becomes nearly a constant) after the initial contact of the two nuclei. A nucleon 
is considered as free if it is not correlated with other nucleon within a spatial distance of 
Ar=3fm and a momentum distance of Ap=300MeV/c. Otherwise, it is bound in a cluster. 
In addition to the nucleon emissions we also calculated all of the fragments during the 
same reaction . Fig. 2 shows the time evolutions of the N„(top windows) and N p (bottom 
windows) for the colliding systems 76 Kr+ 40 Ca, 74 Kr+ 74 Se, 76 Kr+ 76 Kr and 74 Se+ 74 Se with 
neutron-proton ratios 1.04,1.06,1.11 and 1.18 respectively at impact parameter b=4.0 fm 
and beam energy E=150 Me V/ nucleon. To identify the isospin effects of the two-body col- 
lision on the N n (Np) we compare the cases when the isospin dependence of in- medium 
NN cross section is either turned on {a lso ) or off (<j nmso ). Turning on the isospin-dependent 
NN cross section is seen to enhance the momentum dissipation as expected, leading to a 
larger number of proton (neutron) emissions. The isospin effects of the one-body dissipa- 
tion and two-body collsion on the N n (N p ) are identified by using five cases a 
(solid lines), Uo ym + a iso (dot lines), XJ s 2 ym + a iso (dashed lines) ,U{ yrn + a noiso (dashed-dot 
lines), U^™ 1 + a nmso (dashed-dot-dot lines). From Fig. 2 it is clear to see that the varia- 
tions among the values of N p (N ra ) with the same a tso or the same a nmso but different 
symmetry potentials are smaller but the gaps between the N p (N ra ) with a tso and N p (N„) 
with cr nmso are larger, i.e,the values of N p (N n ) depends sensitively on the isospin depen- 
dence of in-medium NN cross section and weakly on the symmetry potentials for the 
neutron-deficient systems. This situation is contrary to the neutron-proton ratio of the 
preeqilibrium emissions in relatively low energy ergions, where the neutron-proton ratio of 
the preeqilibrium emissions depends sensitively on the symmetry potential and weakly on 
the isospin dependence of the in-medium NN cross section 1 19 1. With decreasing beam 
energy the role of two-body collision is reduced. In particular the symmetry potential is 
repulsive for the neutrons and attractive for the protons which tends to make more neu- 
trons than protons unbound i.e. the role of symmetry potentail on the nucleon emissions at 
low energies results in the sizeable differences between neutron emissions and proton emis- 
sions,but the two-body collision produces about the same probability for gaining enough 
energy to become unbound for the the neutrons and protons. As a result,the isospin depen- 
dence of in-medium NN cross section has a small affect on the neutron-proton ratio of the 
preeqilibrium emissions but the role of the symmetry potential on it is strengthened in re- 
gions of relatively low beam energy.But with increasing beam energy the role of two-body 
collisions increases while the role of the mean field is reduced i.e., in the relatively high 
energy region the two-body collision is dominant. Especially in this paper we calculated 
the neutron emissions and the proton emissions not the ratio of them. Finally we can get 
the results as in Fig.2. In Fig.3 is shown the impact parameter-averaged asymptotic values 
of N n (top window) and N p (bottom window) as a function of beam energy for the systems 

76 Kr+ 40 Ca and 76 Kr+ 76 Kr in the fiye 

cases considered as the same as Fig.2. From Fig.3 
it is very clear to see the same conclusion as Fig.2 in the beam energy region from about 
50 to 400 MeV/nucleon. But as the beam energy is decreased to about 50 MeV/nucleon 
the N p (N n ) depends on both the isospin dependence of in-medium NN cross section and 
the symmetry potential.lt is worth mentioning that the minimum energies for remaining 
above property are small different for the systems with the variation of the neutron-proton 
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ratios and mass asymmetries of the colliding systems. 

The above results show strongly that N p (N„) during the neutron-deficient nuclear reac- 
tion can be used as a sensitive probe for extracting information on the isospin dependent 
in-medium NN cross section. Here the above behaviour has been shown for the neutron 
emissions and proton emissions separately,of couse,there is also the same behaviour as a 
probe for the sum of them. 

It is worth mentinoing that in general, the correction effect of the sequential decays on 
the dynamical process of HIC and the number of the nucleon emissions should be consid- 
ered. For example M.B.Tsang et al'studies show that the apparent temperatures measured 
with double ratios of fragment isotope yields display fluctuations that can be attributed 
to the sequential decay of heavier particle unstable nuclei ^ .However the number of nu- 
cleon emissions in our calculations is the final total number of nucleon emissions after 
colliding system has reached at equilibrium which includes the original nucleon emissions 
and sequential decays leading above conclusion. Namely the final total number of nucleon 
emissions dependes sensitively on the isospin effect of in-medium NN cross section and in- 
sensitively on the symmetry potential as a probe for extracting information on the isospin 
dependent in-medium NN cross section in HIC at the relatively high beam energies for 
neutron-deficient colliding systems. Even though the dynamical mechanisms for the origi- 
nal nucleon emissions and the sequential decays are different which produces the correction 
effect on the total number of nucleon emissions but this difference and the correction effect 
of sequential decays on the total number of nucleon emissions do not influence our final 
conclusion. 

4 Summary and conclusions 

In summarywithin the IQMD we studied the isospin effects of one-body dissipation 
and two-body collision on the number of protons (neutrons) emitted during the nuclear 
reaction. The calculated results show strongly that the isospin-dependent in-medium NN 
cross section has a much stronger influence on N p (N n ) but the effects of the symmetry po- 
tential on them are smaller for the neutron-deficent systems in the relatively high energy 
region. Studies of N p (N n ) during the nuclear reactions are proposed to extract informa- 
tion on the isospin-dependent in-medium NN cross section in HIC based on systimatic 
comparisons between the theoretical simulations and experimental studies. 

This work was supported in part by the Major State Research Development Project in 
China under Contract No. G2000077400,I00 person Project of the Chinese Academy of 
Sciences, the National Natural Foundation of China under Grants No.I0I75080,No.l0004012 
and The CAS Knowledge Innovation Project No.KJCX2-SW-N02, 
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Figure captions 

Fig. 1 The correlation between the mean intermediate mass fragment multiplicity N im f 
and the charged particle multiplicity N c . Filled (unfilled) circles represent the exper- 
imental data [36] for the reactions 58 Ni+ 58 Ni( 58 Fe+ 58 Fe) at E=75 MeV/u and the 
solid line (dot line) indicates the IQMD results for 58 Ni + 58 Ni( 58 Fe + 58 Fe).The 
charge number of Aj m j is taken from 3 to 18. 
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Fig. 2 The time evolutions of the N n (top windows) and N p (bottom windows) for the 
colliding systems 76 Kr+ 40 Ca, 74 Kr+ 74 Se, 76 Kr+ 76 Kr and 74 Se+ Se with neutron- 
proton ratios 1.07,1.11,1.11 and 1.18 respectively at impact parameter b=4.0fm and 
beam energy E=150MeV/nucleon for five cases (see text). 

Fig. 3 The impact parameter-averaged values of N n (top window) and N p (bottom win- 
dow)as a function of the beam energies for the systems 76 Kr + 40 Ca and 76 Kr+ 76 Kr 
in the five cases as the same as Fig. 2. 
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